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Abstract 

Iron homeostasis is an important process for flower development and plant fertility. The role of plastids in these processes 
has been shown to be essential. To document the relationships between plastid iron homeostasis and flower biology fur- 
ther, a global study (transcriptome, proteome, metabolome, and hormone analysis) was performed of Arabidopsis flow- 
ers from wild-type and triple atfer1-3-4 ferritin mutant plants grown under iron-sufficient or excess conditions. Some 
major modifications in specific functional categories were consistently observed at these three omic levels, although 
no significant overlaps of specific transcripts and proteins were detected. These modifications concerned redox reac- 
tions and oxidative stress, as well as amino acid and protein catabolism, this latter point being exemplified by an almost 
10-fold increase in urea concentration of atfer1-3-4 flowers from plants grown under iron excess conditions. The mutant 
background caused alterations in Fe-haem redox proteins located in membranes and in hormone-responsive proteins. 
Specific effects of excess Fe in the mutant included further changes in these categories, supporting the idea that the 
mutant is facing a more intense Fe/redox stress than the wild type. The mutation and/or excess Fe had a strong impact 
at the membrane level, as denoted by the changes in the transporter and lipid metabolism categories. In spite of the 
large number of genes and proteins responsive to hormones found to be regulated in this study, changes in the hormo- 
nal balance were restricted to cytokinins, especially in the mutant plants grown under Fe excess conditions. 

Keywords: Ferritin, flowers, functional genomics, iron. 



Introduction 

Flower development and plant fertility are dependent on 
iron homeostasis (Takahashi et al. , 2003). Small organic mol- 
ecules such as nicotianamine (NA) and citrate participate in 



this process by chelating Fe in order to shuttle it through- 
out the plant (Curie et al, 2009; Rellan-Alvarez et al, 2010; 
Schuler et al. , 2012). Overexpressing the barley nicotianamine 
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amino transferase (NAAT) in transgenic tobacco depletes the 
NA pools of these plants, leading to abnormally shaped and 
sterile flowers (Takahashi et al, 2003). NA-Fe(II) complexes 
are transported by members of the YSL gene family (Curie 
et al, 2009; Conte and Walker, 2011). Among them, it is well 
established that the YSL1 and YSL3 genes are required for 
fertility of Arabidopsis flowers, and that pollen function is 
severely reduced in a ysllysB double knockout mutant (Le 
Jean et al, 2005; Waters et al, 2006). Citrate, as an Fe che- 
lator, also appears to be crucial for flower development and 
fertility (Roschzttardtz et al, 2011). The FRD3 gene encodes 
a citrate effluxer (Durrett et al, 2007). Flowers from a null 
frd3 allele have an altered anther structure, and lack pollen 
grain at the anther surface (Roschzttardtz et al, 2011). These 
studies have clearly demonstrated that long-distance Fe traf- 
ficking and its correct allocation within reproductive organs 
are essential for flower development and fertility. 

At the subcellular level, chloroplast Fe homeostasis has 
also recently been linked to flower development and fertility. 
Arabidopsis PIC1 (permease in chloroplastl) is an intrinsic 
membrane protein located at the chloroplast inner envelope, 
and is involved in Fe uptake by plastids (Duy et al, 2007). 
Flowers of PIC1 overexpressor (PIC 1 ox) plants showed 
impaired development, elevated Fe levels, and reprogram- 
ming of their genome expression (Duy et al, 2011). Whereas 
the latter phenotype (PIC 1 ox) resembled that of ferritin 
knockouts (see below), picl knockout mutants accumulated 
ferritin in plastids (Duy et al, 2007). Plant ferritins are Fe 
storage proteins playing a key role in oxidative stress defence 
(Ravet et al, 2009a; Briat et al, 2010a, b). In Arabidopsis, 
flower development of atfer 1-3-4 triple ferritin mutant plants 
is impaired and fertility decreased. Flowers of this mutant 
overaccumulate Fe and have their reactive oxygen species 
(ROS)-scavenging mechanisms activated. Furthermore, the 
absence of ferritin in this mutant led to a strong deregulation 
of the expression of 24 genes encoding metal transporters 
in flowers and stems, including the PIC1 gene (Ravet et al. , 
2009a). 

The extent to which the impact of Fe homeostasis on plant 
development is communicated by plant hormones also needs 
to be considered. Indeed, the interactions between Fe metabo- 
lism and plant hormones have received increased attention in 
the last few years (Romera et al, 2011). Auxin, ethylene, and 
nitric oxide (NO) increased in Fe-deficient plants in relation 
to the up-regulation of root Fe acquisition genes (Bacaicoa 
et al, 2011; Lingam et al, 2011), whereas jasmonate was 
found to suppress this gene induction (Maurer et al, 2011). 
Auxin also plays a role in lateral root elongation regulated by 
local Fe concentrations (Giehl et al, 2012), and abscisic acid 
(ABA) plays a role in the regulation of some ferritin genes in 
response to Fe excess (Lobreaux et al, 1993). 

The main aim of this study was to document the roles that 
Fe homeostasis plays in flower development and fertility, in 
particular when plastid Fe homeostasis is altered. To attain 
this goal, transcriptomic, proteomic, and metabolomic analy- 
sis were performed, and hormone concentrations were deter- 
mined, on the same samples of Arabidopsis flowers collected 
from wild-type (WT) and triple atferl-3-4 mutant plants 



grown under standard or Fe excess conditions. This approach 
enables the description of the impact of the perturbed Fe 
homeostasis caused by the loss of ferritin and/or iron over- 
load at various molecular levels in flowers. In addition, the 
study provides a start for generating a coherent scheme of the 
role that ferritins play in flower biology. 

Materials and methods 

Plant materials, growth conditions, and sampling of flowers 

Arabidopsis Columbia (Col-0) and the previously described mutant 
line atferl-3-4 (Ravet et al, 2009a) were sown on soil (Humin sub- 
strate N2) in pots as single plants and kept at 4 °C for 2 d prior to 
being grown in a growth chamber under long-day conditions (16h 
light/8 h dark, 23 °C). Plants were irrigated once a week with reverse 
osmosis water without Fe added (standard soil) (WT and atferl-3- 
4), or supplemented with Fe excess (2mM Fe-EDDHA) (WT++ and 
atferl-3-4++). 

Flowers were harvested as a mix of development stages 1 and 2, 
as previously defined (Ravet et al, 20096), six times every 2 d on 
10 individuals from 4- week-old atferl-3-4 mutant and WT plants. 
Flowers were immediately frozen in liquid nitrogen after harvest. 
To guarantee comparability of data, all flower samples were pro- 
duced under the same conditions in the same growth chamber of 
only one laboratory (Montpellier). Frozen flower samples, RNA, 
and proteins were then dispatched from the Montpellier laboratory 
to other laboratories (Miinchen, Zaragoza, Davis, and CIPAV) for 
transcriptomics, proteomics, metabolomics, and hormone dosage, 
respectively. 



RNA extraction and transcriptome analysis 

Flowers from different plants (collected three times from 10 indi- 
viduals each) and grown with and without Fe excess were ground 
in liquid nitrogen. RNAs were extracted using an RNeasy Mini 
Kit (Qiagen). For microarray analysis, 5 ug of RNA from flowers 
of the atferl-3-4 mutant and the WT were processed and hybrid- 
ized to Affymetrix GeneChip Arabidopsis ATH1 Genome Arrays 
as described in Duy et al. (2011). For quality check and normali- 
zation, raw intensity values were processed with Robin software 
(Lohse et al, 2010) default settings. Specifically, for background 
correction, the robust multiarray average normalization method 
(Irizarry et al., 2003) was performed across all arrays (between- 
array method). Statistical analysis of differential gene expression of 
atferl-3-4 versus WT flowers was carried out using the linear model- 
based approach developed by Smyth (2004). Thereby, differential 
clustering of n=2 replicates, each representing pooled tissue of 10 
individual plants, was obtained. The obtained P-values were cor- 
rected for multiple testing using the nestedF procedure, applying a 
significance threshold of 0. 10 for mutant versus the WT and 0.20 for 
Fe excess versus Fe-sufficient comparisons, in combination with the 
Benjamini and Hochberg (1995) false discovery rate control. To ver- 
ify differential expression, the transcript content of 12 representative 
genes, classified as regulated by microarrays (Tables 1-3), was re- 
analysed by quantitative real-time reverse transcription-PCR (qRT- 
PCR; Supplementary Table SI available at JXB online). Because 
here only one gene for one comparison showed deviating behaviour, 
micorarray data were considered as robust. 



Protein extraction 

Flowers from different plants were ground in a mortar in liquid 
nitrogen, and the ground material was stored at -80 °C. The flower 
powder (~1 g) was resuspended in 2 ml Eppendorf tubes using 1 ml 
of 10% trichloroacetic acid (TCA), 0.07% p-mercaptoethanol in 
acetone, incubated for 45min at -20°C (shaking at 5, 10, and 
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15min), and centrifuged for 5min at 18 000 rpm at 4 °C. The 
supernatant was discarded and the pellet was resuspended in 
1.5 ml of 0.07% p-mercaptoethanol in acetone, vortexed, incu- 
bated for 5min in an ultrasonic bath, vortexed again, cooled 
with liquid N 2 , incubated for 30min at -20 °C, and centrifuged 
for 5min at 18 000 rpm at 4 °C. The supernatant was discarded 
and the pellet was dried in a SpeedVac for lOmin and stored at 
-20 °C until analysis. 

Pellets were solubilized in 250 pi of sample rehydration buffer 
containing 8M urea, 2% (w/v) CHAPS, 50 mM dithiothreitol 
(DTT), 2mM phenylmethylsulphonyl fluoride (PMSF), and 0.2% 
(v/v) IPG buffer pH 3-10 (GE Healthcare, Uppsala, Sweden), 
stirred at 1000 rpm for 3 h at 38 °C in a Thermomixer Comfort device 
(Eppendorf AG, Hamburg, Germany), centrifuged at 10 000 g for 
15min at room tempertare, and filtered with 0.45 um ultrafree-MC 
filters (Millipore, Bedford, MA, USA). Protein concentration was 
quantified immediately after extraction using Bradford's method 
with bovine serum albumin (BSA) as standard (Sigma). 

In the proteomic and metabolomic analysis, a given genotype per 
treatment combination will be referred to hereafter as a class. 



Two-dimensional electrophoresis (2-DE) analysis 

Protein profiles were obtained using five independent biological rep- 
licates per treatment. The first dimension immunoelectric focusing 
(IEF) separation was carried out on 7cm ReadyStrip IPG Strips 
(BioRad) with a 4-7 linear pH gradient in a Protean IEF Cell 
(BioRad), and the second dimension SDS-PAGE in 8 x 10 x 0.1 cm 
gels, as described in detail elsewhere (Rodriguez-Celma et al., 
2011). Gels were subsequently stained with Coomassie blue G-250 
(Serva, Barcelona, Spain). 2-DE gel image analysis was carried out 
as described in Rodriguez-Celma et al. (2013). Briefly, stained gels 
were scanned, and experimental M r and pi values calculated by 
mobility comparisons with standard markers. Spot detection, gel 
matching, and analysis were performed with PDQuest 8.0 software 
(BioRad). Missing spot volumes were estimated using a sequen- 
tial K-Nearest Neighbor algorithm. Consistent spots were those 
present in at least 80% of the replicates from at least one class. 
Differentially expressed spots were defined using a Student's ?-test 
value of P < 0.10. Protein response ratios were defined as the abun- 
dance in a given class divided by the abundance in the control class; 
when ratios were <1, the inverse was taken and the sign changed. 
Only proteins with mean response ratios >2.0 or below -2.0 were 
considered relevant and analysed further. Partial least square (PLS) 
analysis was applied to the identified protein species using SPSS 
v. 15 (SPSS Inc., Chicago, IL, USA). 



Protein in-gel digestion and identification by nano-high 
performance liquid chromatography-tandem mass 
spectrometry 

Selected spots were excised, in-gel digested with trypsin, and analysed 
by nano-HPLC-tandem mass spectrometry as described in detail in 
Rodriguez-Celma et al (2013). Protein identification was performed 
by searching in the non-redundant NCBI 20120324 (17 612 906 
sequences; 6 046 335 507 residues) and Plants_EST EST_111 (155 120 
040 sequences; 27 373 983 152 residues) databases using the Mascot 
search engine (Matrix Science, London, UK). The search parameters 
were: monoisotopic mass accuracy, peptide mass tolerance ±0.2 Da, 
fragment mass tolerance ±0.6Da, one allowed missed cleavage, 
allowed fixed modification carbamidomethylation (Cys), and variable 
modification oxidation (Met). Positive identification was assigned with 
Mascot scores above the threshold level (P < 0.05), at least two iden- 
tified peptides with a score above homology, 10% sequence coverage, 
and similar experimental and theoretical molecular weight and pi. The 
GO biological process annotation (http://www.geneontology.org) of 
the individual identified proteins was used for classification. 
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Metabolite analysis and data processing 

Metabolite profiles were obtained using six independent biological 
replicates per treatment. Metabolite extraction was carried out as 
previously described for leaf extracts (Fiehn and Weckwerth, 2003; 
Fiehn et al, 2008; Rellan- Alvarez et al, 201 1). Dried extracts were 
derivatized as described elsewhere (Fiehn et al. , 2008). Derivatized 
samples (1 ul) were injected randomly in split-less mode with a cold 
injection system (Gerstel, Mulheim an der Ruhr, Germany) and 
separated in a GC device (Agilent 6890, San Jose, CA, USA) using 
an integrated guard column (Restek, Bellefonte, PA, USA) and an 
Rtx 5Sil MS column (30 mx0.25mm, 0.25 urn film thickness). The 
GC device was connected to a Leco Pegasus IV time-of-flight mass 
spectrometer (TOFMS) controlled with Leco ChromaTOF software 
v.2.32 (Leco, St. Joseph, MI, USA). Peak detection and mass spectra 
deconvolution were performed with Leco Chroma-TOF software 
v.2.25, and GC-MS chromatograms were processed as described 
previously (Fiehn et al, 2008). 

Metabolite data were normalized using the sample fresh weight and 
the sum of all metabolite peak heights in a single run, to account for 
small GC injection variations. The resulting data were multiplied by 
a constant factor in order to obtain values without decimal figures. 
Peaks were identified using the Fiehn Lib database (http://fiehnlab. 
ucdavis. ed u/ proj ects/FiehnLib/index_html) . Metabolite spectra, 
including those of identified metabolites and unknowns, can be que- 
ried and downloaded at the Fe Chlorosis Database (fechlodb.fiehnlab. 
ucdavis.edu). Data were analysed to check for possible correlations 
between peak height values and peak variance, and since a positive 
correlation was found a log 10 transformation of the data was carried 
out to avoid variance-mean dependence (Chich et al. , 2007). 

Statistical analysis of the normalized data was carried out with 
Statistica software (v.9.0. StatSoft, Inc., Tulsa, OK, USA). Only 
those metabolites present in at least 80% of the samples from at least 
one class were considered. Significant changes in metabolite levels 
were detected for each class using one-factor analysis of variance 
(ANOVA; P < 0.05). Metabolite response ratios were defined as the 
level in the class divided by the level in the control class; when ratios 
were <1, the inverse was taken and the sign changed. Only metab- 
olites with mean response ratios (using the non-log 10 transformed 
data) >1.5 or below -1.5 were considered relevant and are discussed 
in this study. Multivariate analysis (supervised PLS) was used to 
study the clustering of the samples, as well as to find out which 
metabolites were responsible for the separation between classes. 

For all the metabolomic profiling analyses, the recommended 
guidelines provided by Fiehn et al. (2008) were followed. Full 
documentation of metabolite profiling data acquisition and inter- 
pretation, according to Fernie et al. (2011), is also disclosed in 
Supplementary Table S2 at JXB online. 



Hormone analysis 

The concentration of the main plant hormones in the flower sam- 
ples of WT and atferl-3-4 Arabidopsis plants was determined using 
HPLC-electrospray-mass spectrometry (HPLC-ESI-MS/MS) as 
described below. 

The following hormones were analysed and quantified: indole- 
3-acetic acid (IAA), ABA, trans-zeatm (Z), trans- and cw-zeatin 
riboside (t-ZR and c-ZR), dihydrozeatin riboside (DHZR), and 
isopentenyladenosine (iPR). 

The extraction and purification of IAA and ABA were carried out 
using the method described in detail by Bacaicoa et al. (201 1), using 
250 mg of frozen flowers and 4 ml of extraction medium. Overnight 
extraction, centrifugation, re-extraction, and purification of the 
extract were done as in Bacaicoa et al. (201 1). Before the injection in 
the HPLC-ESI-MS/MS system, the solution was centrifuged at 12 
000 g for 5 min. IAA and ABA were determined by HPLC-ESI-MS/ 
MS using multiple reaction monitoring (MRM) in the negative-ion 
mode, employing multilevel calibration curves with deuterated hor- 
mones as internal standards, as indicated in Bacaicoa et al. (2011). 



Compound-dependent parameters are listed in Supplementary 
Table S3 at JXB online. The extraction and purification of cyto- 
kines were carried out using the method described by Dobrev and 
Kaminek (2002), with some modifications: 250 mg of frozen flower 
tissue were ground in a mortar with liquid N 2 and homogenized 
with 4ml of pre-cooled (-20 °C) methanol: watenformic acid (15:4:1, 
v/v/v). Deuterium-labelled cytokinin internal standards ([ 2 H 5 ]Z, 
[ 2 H 5 ]tZR, [ 2 H 6 ]iP, [ 2 H 6 ]iPR, [ 2 H 7 ] BA, [ 2 H 7 ]BAR, and [ 13 C 5 ]oT; all 
from Olchemim Ltd, Olomouc, Czech Republic), were added (40 ul 
of a stock solution of 50 ng ml -1 of each standard in methanol) to 
the extraction medium. Extraction was then performed according 
to Bacaicoa et al. (2011) until reaching the aqueous phase after 
evaporation at 40 °C. Then, 2 ml of 1M formic acid was added, 
and applied to an Oasis MCX column (3cm 3 , 60mg) (Waters Co., 
Milford, MA, USA; Ref. 186000254), pre-conditioned with 4ml of 
methanol and 2 ml of 1 M formic acid. The column was washed suc- 
cessively with 2 ml of 1 M formic acid, 2 ml of methanol, and 2 ml of 
0.35 M NH 4 OH, and the cytokinin bases, ribosides, and glucosides 
were eluted with 2ml of 0.35 M NH 4 OH in 60% (v/v) methanol. This 
eluted fraction was evaporated to dryness in a vortex evaporator and 
re-dissolved in 250 ul of methanol:0.05% formic acid (40:60, v/v). 
Before the injection in the HPLC-ESI-MS/MS system, the solution 
was centrifuged at 12 000 g for 5min. Cytokinins were quantified by 
HPLC-ESI-MS/MS as in Bacaicoa et al. (201 1), but with a different 
reverse-phase column (Tracer Excel 120 ODSA 3 um, 100 x 4.6mm; 
Teknokroma, Barcelona, Spain). A linear gradient of methanol 
(A) and 0.05% formic acid in water (B) was used: 35% to 95% A in 
1 1 min, 95% A for 3 min, and 95% to 35% A in 1 min, followed by a 
stabilization time of 5min. The flow rate was 0.25ml mhr 1 , the injec- 
tion volume was 40 ul, and column and sample temperatures were 
30 °C and 20 °C, respectively. The detection and quantification of 
cytokinins were carried out using multiple MRM in the positive-ion 
mode, employing multilevel calibration curves with deuterated cyto- 
kinins as internal standards. Compound-dependent parameters are 
listed in Supplementary Table S4 at JXB online. The source param- 
eters were: curtain gas, 172.37 kPa; GS1, 344.74 kPa; GS2, 413.69 
kPa; ion spray voltage, 5000 V; and temperature, 600 °C Data sam- 
ples were processed using Analyst 1.4.2 Software from AB SCIEX 
(Framingham, MA, USA). 

Results 

Comparison of the flower transcriptomes of Col-0 
and atferl -3-4 Arabidopsis plants grown under iron- 
sufficient or excess conditions 

Differential gene expression in flowers of the triple ferritin 
mutant atferl-3-4 (Ravet et al, 2009a) and Col-0 WT plants, 
grown in Fe-sufficient control conditions or under excess Fe 
(atferl-3-4++ and WT++, respectively), was monitored by 
microarray analysis. The differential regulation observed was 
low, with a total of 276 genes significantly regulated (Fig. 1A; 
Supplementary Table S5 at JXB online). Most fold changes 
ranged between 1.4 and 3.5, and only a few showed strong 
alterations (7). All ferritin genes mutated in atferl-3-4 (Ravet 
et al. , 2009a) were down-regulated when compared with the 
WT, with strong effects for the knockouts of AtFERl and 
AtFER4 and a mild decrease for the knockdown of AtFER3 
(Fig. 1C), the latter only being significant when the atferl- 
3-4 mutant was challenged with excess Fe (see number 3 in 
Supplementary Table S5). When flower transcript levels in 
mutants versus the WT in all three possible combinations were 
compared, 121 {atferl-3-4 versus WT), 187 (atferl-3-4++ 
versus WT++), and 146 (atferl-3-4++ versus WT) genes were 
found to be significantly regulated (adjusted P-value <0.10; 
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(A) atfer1-3-4 vs. WT atfer1-3-4++ vs. WT++ (B) atfer1-3-4++vs. atfer1-3-4 WT++ vs. WT 
121:^95^26 187:{}88f99 28: {} 6 J 22 25:{}l4^11 




atfer1-3-4++ vs. WT 

„^ A, r l„ FER4 89.6 91.3 1606.3 1413.1 

146: if 75 T 71 

Fig. 1. Changes in Arabidopsis flower transcripts of atferl -3-4 mutants and the wild type. Overlapping regulation is shown by Venn 
diagrams. Numbers of regulated genes are given; up- and down-regulation is indicated by orange and green arrows, respectively 
(for fold changes see Supplementary Table S5 at JXB online). (A) Differentially expressed genes (P-value <0.1) in flowers of atferl '-3-4 
mutants compared with the Col-0 wild type (WT) under Fe-sufficient and Fe excess (++) conditions. In total, 276 genes were significantly 
regulated when comparing atferl -3-4 versus WT, atferl-3-4++ versus WT++, and atferl-3-4++ versus WT. Gene groups of interest 
discussed in the text are highlighted by a grey background. (B) Transcript changes (P-value <0.2) induced by Fe excess (++) in atferl -3- 
4 and the wild type using the respective Fe-sufficient conditions as control. Note that all genes regulated in atfer1-3-4++ versus atferl- 
3-4 and WT++ versus WT are also differentially expressed in one of the comparisons depicted in (A). (C) Mean signal values (arbitrary 
units) of mRNA from mutated ferritin genes (AtFERl , AtFER3, and AtFER4) in flowers of all samples analysed. 



Fig. 1A). Interestingly, in control conditions, the majority of 
genes (95 of 121) in atferl-3-4 were up-regulated, whereas 
upon Fe excess, the increase and decrease of transcript num- 
bers was more evenly distributed, when compared both with 
WT++ (88 up and 99 down) and with WT (75 up and 71 
down). When compared with transcript regulation in flowers 
overexpressing the chloroplast Fe permease PIC1 (PIC 1 ox; 
Duy et al, 2011), a high number of overlapping genes was 
found (63 out of 276, 23%; see highlighted genes in Tables 1-3 
and Supplementary Table S5). In comparison with mutant 
versus WT changes, regulation in response to Fe was less pro- 
nounced and less significant (adjusted P-value <0.20), with 
only 28 (atferl-3-4++ versus atferl-3-4 ) and 25 (WT++ ver- 
sus WT) genes showing differences (Fig. IB). Iron excess in 
flowers of the atferl-3-4 mutant mainly induced down-regu- 
lation of genes (22 of 28), while in the WT the increases and 
decreases in transcript numbers were evenly distributed ( 14 up 
and 1 1 down). To correlate data output from transcriptomics, 
proteomics, and metabolomics throughout the entire study, 
the following 17 functional categories were chosen for sort- 
ing of differentially expressed genes and proteins (see Tables 
1-4 and Supplementary Table S5): metal (iron) homeostasis 
and binding, amino acid metabolism, protein metabolism 
and modification, carbohydrate metabolism, photosynthe- 
sis, cofactor biosynthesis, lipid metabolism, stress and redox, 
DNA/RNA processes, cell cycle and cell biogenesis, cell wall 
modification, secondary metabolism, hormone metabolism, 
other processes, signalling, transport, and unknown function. 
For all putative functional annotations cited in the following 



sections: (i) the MAPMAN Arabidopsis pathway analysis 
program (Thimm et al, 2004); (ii) TAIR (The Arabidopsis 
Information Resource, TAIR Gene Model on www.arabidop- 
sis.org, 16 June 2012); and (hi) the ARAMEMNON database 
(aramemnon.botanik.uni-koeln.de; Schwacke et al, 2003) 
were referred to. 

The focus was first on overlaps in mutant versus WT sam- 
ples in all three comparisons {atferl-3-4 versus WT, atferl- 
3-4++ versus WT++, and atferl-3-4++ versus WT). This 
overlap contained 54 differentially expressed genes (44 with 
P-value <0.05), which reflects the influence of the ferritin null 
mutant background regardless of the Fe status of the plants 
(Fig. 1 A). Disposition of these genes into functional categories 
(Table 1) revealed that the biggest focus of transcript regula- 
tion was on stress- and redox-related functions, accounting for 
28% of the total changes. In this category, all genes were up- 
regulated in atferl-3-4, and most of them (7 out of 15) were 
Fe-haem-binding oxidoreductases. Further, besides the 
AtFERl and AtFER4 genes, the largest fold changes (>10.0) 
were also found in this category: an FAD-binding berberine 
family protein (Atlg26390), and, when plants were challenged 
with excess Fe, a P450 monooxygenase (At4g31970). In addi- 
tion to stress and redox, transcriptional regulation and hor- 
mone metabolism were additional targets for gene regulation 
(i.e. more than three genes differentially expressed). In gen- 
eral, hormone-related genes were up-regulated, whereas the 
transcriptional regulation category contained genes up- and 
down-regulated. For only three genes, Fe excess treatment in 
the comparisons atferl-3-4++ versus WT and atferl-3-4++ 
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versus WT++ induced a reciprocal regulation (i.e. transcript 
decrease) when compared with Fe-sufficient conditions in 
atferl-3-4 versus WT (increase of transcript content; Table 1): 
TCH4 (At5g57560), a xyloglucan endotransglycosylase 
involved in cell wall modification; At4g27280, a Ca-binding 
EF hand family protein; and At4g36500, a protein of 
unknown function predicted to be localized in mitochondria 
and regulated in response to various stresses or auxin (TAIR 
annotation). The largest overlap with gene regulation in flow- 
ers overexpressing the chloroplast Fe permease PIC1 (PIC 1 ox; 
Duy et al, 2011) was for genes involved in hormone metabo- 
lism. To map the effect of excess Fe on the mutant versus WT 
regulation, those genes that were differentially expressed in 
the atferl-3-4++ versus WT++ comparison were selected (74 
genes, 14 with a P-value <0.05; Fig. 1A, Table 2). In addition 
to the functional categories described above, protein metabo- 
lism and modification, cell wall modification, and transport 
were also extra targets for gene regulation. Metal homeosta- 
sis-associated genes were represented by the ferric reductase 
oxidase FR03 (Atlg23020), induced by Fe deficiency, and 
involved in Fe transport processes in the root vascular cylinder 
(Mukherjee et al, 2006), which was increased in the mutant 
when compared with the WT. Overall, genes related to protein 
and hormone metabolism were down-regulated, while those 
belonging to stress and redox processes, in contrast to those 
listed in Table 1 (overlaps of all mutant- WT comparisons), 
were now evenly distributed among increases and decreases. 
With regard to the new targets, protein metabolism and cell 
wall modification showed a general down-regulation, whereas 
transport-related genes were evenly distributed (three up- and 
three down-regulated). Among the transporters increased in 
atferl-3-4++ when compared with WT++, two related to Al 
tolerance were found: the phloem-localized ABC transporter 
permease-like protein ALS3 ( At2g37330) (Larsen et al. , 2005) 
and ALMT1 (Atlg08430), an Al-activated malate trans- 
porter (Liu et al, 2009, 2012). The other up-regulated gene 
was CAT6 (At5g04770), a putative cationic amino acid trans- 
porter. Down-regulated transporter genes were: the ABC- 
transporter ABCG10 (Atlg53270); DTX41 (At3g59030), 
a putative vacuolar flavonoid/H+ antiporter of the MATE 
family; and a Secl4p-like phosphatidylinositol transfer fam- 
ily protein (Atlg30690) with potential a-tocopherol trans- 
port function. The 43 genes differentially expressed only in 
atferl-3-4++ versus WT (14 with a P-value <0.05; Fig. 1A, 
Table 3) should mainly reflect the effect of Fe excess on the 
mutant only. The targets for gene regulation within this com- 
parison were again stress and redox, DNA/RNA processes, 
and hormone metabolism, and in general showed a down- 
regulation. However, the largest change was an up-regulation 
(8-fold increase) of ACCD (AtCg00500), the plastid-encoded 
carboxytransferase |3-subunit of the acetyl-CoA carboxylase 
complex. This complex catalyses the carboxylation of acetyl- 
CoA to produce malonyl-CoA, the first committed step in 
fatty acid synthesis. Chloroplast-encoded genes appeared to 
be predominantly up-regulated in the atferl-3-4++ versus WT 
comparison. Besides ACCD, YCF5, which is involved in the 
assembly of the cytochrome c complex (AtCg01040, Table 3, 
4-fold increase), and RPS16, a plastid-localized S16 ribosomal 



protein (AtCg00050, Table 1, 14-fold increase) were also found 
to be strongly induced in atferl-3-4++ versus WT. This sug- 
gests that under excess Fe the loss of ferritin function affects 
chloroplast gene expression and in consequence metabolism. 

Comparison of the flower proteomes of Col-0 and 
atferl -3-4 Arabidopsis plants grown under iron- 
sufficient or excess conditions 

Typical real scans of 2-DE gels obtained from protein 
extracts of flowers from Col-0 and triple ferritin mutant 
grown in control conditions or under Fe excess (WT, WT++, 
atferl-3-4, and atferl -3-4++, respectively) are shown in 
Fig. 2A-D. The average number of detected spots (mean 
±SD; n=5) was 362±35, 328±39, 341 ±25, and 350±40 
in WT, WT++, atferl-3-4, and atferl-3-4++, respectively 
(Supplementary Table S6 at JXB online). The total number 
of consistent spots was 377. A composite averaged virtual 
map containing all spots present in all classes is shown in 
Fig. 2E. 

A first comparison was made taking the WT protein profile 
as control. Approximately 3% (11), 6% (20), and 6% (21) of 
the spots in the WT++, atferl-3-4, and atferl-3-4++ classes, 
respectively, showed significant relative abundance changes 
(spots are highlighted in Fig. 2F-H; Supplementary Table S6 
at JXB online). When a second comparison was carried out 
with the same criteria between atferl-3-4++ and atferl-3-4, 
6% (22) of the spots showed changes (Fig. 3A; Supplementary 
Table S6). A third comparison between atferl-3-4++ and 
WT++ revealed that 3% (12) of the spots showed changes 
(Fig. 3B; Supplementary Table S6). The total number of spots 
having consistent changes between classes was 63, and 98% 
of them (62) were unambiguously identified by MS (Table 4; 
Supplementary Fig. S1A-E). The PLS analysis showed a 
good separation between treatments WT and atferl-3-4 ++, 
whereas WT++ and atferl-3-4 clustered together and apart 
from the other two (Supplementary Fig. S2A). 

The Fe excess treatment in WT plants only caused 
decreases in the relative signal intensity of 11 spots (green 
symbols in Fig. 2F) and all of them were identified (Table 4; 
Supplementary Fig. S1A at JXB online). Proteins belonged 
to different metabolic pathways, including amino acid and 
protein metabolism (spots 3, 6, and 13), carbohydrate metab- 
olism (spot 25), photosynthesis (spot 28), cofactor biosynthe- 
sis (spot 34), stress-related processes (spots 39 and 40), DNA 
replication (spot 43), and N metabolism (spot 53). Spot 56 
was an uncharacterized protein. The largest fold decreases 
in this comparison were found for an elongation factor 
involved in protein translation (25-fold; spot 13) and for a 
(3-ureidoropionase (10-fold; spot 53), a hydrolase acting in 
pyrimidine catabolism that yields (3-alanine, C0 2 , and NH 3 . 

When comparing the averaged map of flowers from atferl- 
3-4 versus WT, two spots increased in signal intensity (orange 
symbols in Fig. 2G). They were identified as the large subunit 
of Rubisco and an uncharacterized protein (spot 27 and spot 
58, respectively, in Table 4 and Supplementary Fig. SIB at 
JXB online). Seventeen spots decreased, whereas one more 
(spot 50) was no longer detected (green and blue symbols, 
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Fig. 2. 2-DE IEF-SDS PAGE protein profile maps of whole flower extracts from Arabidopsis thaliana plants, and changes in spot 
abundance taking the WT as control. Proteins were separated in the first dimension in linear IPG pH 4-7 gel strips and in the second 
dimension in 1 2% acrylamide vertical gels. Scans of typical gels of flower extracts from WT, WT++, atfer1-3-4, and atfer1-3-4++ plants 
are shown in A, B, C, and D, respectively. To facilitate visualization of the studied spots, a virtual composite image (E) was created 
containing all spots present in the real gels A, B, C, and D. F, G, and H show spots whose relative abundance changed in WT++, atferl- 
3-4, and atferl-3-4++ when compared with the WT. Spots that decreased in abundance or were no longer detected when compared 
with control maps are marked with green and blue symbols, respectively, and those with increased abundance or newly detected spots 
are marked with orange and red symbols, respectively. 



respectively, in Fig. 2G). All of them were identified (Table 4; 
Supplementary Fig. SIB) and assigned to amino acid (spots 
4 and 5), carbohydrate (spots 22-24 and 26), and protein 
(spots 15 and 16) metabolic processes, cofactor biosynthesis 
(spot 34), lipid metabolism (spot 36), stress (spots 38 and 41), 
DNA- and RNA-related processes (spots 42-44), cell cycle 
(spot 47), and secondary metabolism (spot 50). Spot 59 was 
an uncharacterized protein. The largest decreases found were 
for the GSL lipase, an extracellular enzyme involved in lipid 
degradation (24-fold, spot 36), the monothiol glutaredoxin- 
S15 (15-fold, spot 38), a mitochondrial enzyme containing 
an Fe-S cluster involved in cell redox homeostasis, and a 
strictosidine synthase (no longer detected, spot 50), which is 
involved in indole alkaloid biosynthesis. 

Comparison of atferl-3-4++ versus WT showed eight 
spots that increased and one appearing de novo (orange and 
red symbols in Fig. 2H, respectively). All of them were iden- 
tified (Table 4; Supplementary Fig. SIC at JXB online) and 
five were ascribed to amino acid and protein (spots 1, 10, and 
11) and carbohydrate metabolism (spots 18 and 21). Spot 52 
was in the Others category, and spots 57, 60, and 61 were 
uncharacterized proteins. Ten spots decreased (green symbols 
in Fig. 2H), and two were no longer detected (blue symbols in 
Fig. 2H). All of them were identified (Table 4; Supplementary 
Fig. SIC) and assigned to protein metabolism (spots 13, 14, 
and 17), photosynthesis (spots 29-31), lipid metabolism (spot 
36), stress (spot 38), RNA-associated processes (spot 42), 
cell cycle (spot 46), and others (spots 54 and 55). The larg- 
est changes within this comparison were found for an elonga- 
tion factor (spot 13), an RNA-binding protein related to the 



stress response (spot 42), which was no longer detected, and 
the cytosolic fructose- 1,6-biphosphatase (spot 18), which was 
detected de novo. 

In the atferl-3-4++ versus atferl-3-4 comparison, 18 spots 
increased and two appeared de novo (orange and red sym- 
bols in Fig. 3A, respectively). Nineteen of them were identi- 
fied (Table 4; Supplementary Fig. SID at JXB online) and 
ascribed to amino acid (spots 1 and 5), protein (spots 7, 9 
and 12), lipid (spot 35), carbohydrate (spots 18-20 and 26), 
and secondary (spot 50) metabolism, thiamine biosynthesis 
(spots 32 and 33), stress (spot 41), and cell cycle (spots 45 and 
47). Spot 52 was in the Others category, and spots 61 and 62 
were uncharacterized proteins. Spot 63 could not be identi- 
fied. One spot, 49, was no longer detected and was identified 
as (3-xylosidase 4, a protein assigned to cell wall modifica- 
tion (green and blue symbols in Fig. 3 A; Table 4; Fig. SID), 
and another spot, 56, decreased in relative intensity and was 
an uncharacterized protein. In addition to the proteasome 
subunit (3-type-2-B (spot 7) and the strictosidine synthase 
(spot 50), which were detected de novo, the largest change in 
this comparison was an increase in a chloroplastic haloacid 
dehalogenase-like protein (21-fold, spot 62). 

The atferl-3-4++ versus WT++ comparison revealed 
10 spots increasing in intensity, and one appearing de novo 
(orange and red symbols in Fig. 3B, respectively). All of 
these were identified (Table 4; Supplementary Fig. S1E at 
JXB online) and included in amino acid (spots 1-3), protein 
(spot 8), carbohydrate (spots 18 and 20), and secondary (spot 
51) metabolism, stress (spot 40), and cell wall (spot 48) cat- 
egories. Spot 52 was in the Others category, and spot 61 was 
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Fig. 3. Changes in spot abundance in the 2-DE IEF-SDS-PAGE protein profile maps of whole flower extracts from after? -3-4++, taking 
atferl-3-4 as the control. Spots that decreased or were no longer detected when compared with control maps are marked with green 
and blue symbols, respectively, and those with increased response ratios (i.e. abundance) or newly detected spots are marked with 
orange and red symbols, respectively. 



an uncharacterized protein. One decreasing spot, identified 
as a GDSL esterase/lipase, was related to lipid metabolism 
(green symbol in Fig. 3B; spot 37 in Supplementary Fig. 
S1E). The largest changes within this comparison were a 
100-fold increase in a cysteine protease (spot 8) and a 10-fold 
increase in a xyloglucan endotransglucosylase (spot 48) that 
participates in cell wall construction of growing tissues by 
modifying xyloglucan polymers, an essential constituent of 
the primary cell wall. 

The information on the changes in protein abundance 
overlapping in the different classes was demonstrated using 
Venn diagrams, as for transcriptomic data analysis (Fig. 4). 
No protein overlap was found in the atferl-3-4 versus WT, 
atferl-3-4++ versus WT++, and atferl-3-4++ versus WT 
comparisons. The second group of interest, containing pro- 
teins changing only in the atferl-3-4++ versus WT++ com- 
parison (Fig. 4, Table 4), contained eight protein species, 
most of them (seven) showing increases in relative abun- 
dance, except for one (spot 37, identified as a GDSL esterase/ 
lipase). The amino acid/protein metabolism category, which 
contained a D-3-phosphoglycerate dehydrogenase, an argin- 
ase, and a cysteine protease (spots 2, 3, and 8, respectively), 
was a target for regulation at the protein level. In agreement 
with the target categories found in the transcriptomic analysis 
for gene regulation, protein species in this group were related 
to stress and redox, including an Fe oxidoreductase (flavonol 
synthase, spot 51), a galactose oxidase (spot 20), and a glu- 
tathione 5-transferase (spot 40), and to cell wall modification 
(xyloglucan endotransglucosylase, spot 48). Although the 
hormone category did not contain any protein species, at least 
three of the proteins in this second group (spots 40, 48, and 
51) are auxin regulated (UniProt DB annotation). The third 
group of interest contained 14 protein species changing only 
in the atferl-3-4++ versus WT comparison (Fig. 4, Table 4). 
This analysis revealed that protein metabolism (spots 10, 11, 
13, 14, and 17) and photosynthesis [two spots for the cytb6f 
protein (spots 29 and 30) and OEC23 (spot 31)] were target 
categories for regulation at the protein level. 



Comparison of the flower metabolomes of Col-0 
and atferl -3-4 Arabidopsis plants grown under iron- 
sufficient or excess conditions 

The metabolite profile of whole flower extracts from the 
A. thaliana triple ferritin mutant atferl-3-4 and the Col-0 WT, 
grown in control conditions and under excess Fe (atferl -3- 
4++ and WT++, respectively), were studied using GC-MS. 
Only metabolites present in at least 80% of any sample class 
(in five out of six replicates) were taken into consideration. 
A total of 251 of such consistently present metabolites were 
detected, and 103 of them were identified using the Fiehn Lib 
database, whereas 148 were unknowns (a complete metabolite 
list is shown in Supplementary Table S7 at JXB online). 

When comparing the metabolite profiles of the different 
sample classes with those of the WT, 10 identified metabo- 
lites showed significant changes (at P < 0.05) and mean 
response ratios >1.5 or below -1.5. One, three, and 10 of the 
known metabolites showed changes in the WT++, atferl-3- 
4, and atferl-3-4++ classes, respectively (Table 5). The larg- 
est increases in abundance were for urea (10-fold increase in 
atferl -3-4++), 2-hydroxyglutaric and threonic acids (2- to 
3-fold increases in atferl-3-4 and atferl -3-4++), and gluconic 
acid lactone (2-fold increases in atferl -3-4++). The largest 
decreases in abundance were for ascorbic acid (80% decreases 
in WT++ and atferl -3 -4++), and salicin (75% in atferl-3- 
4++). On the other hand, significant changes were observed 
in a total of 25 unknown metabolites: three, seven, and 22 
of them in the WT++, atferl-3-4, and atferl-3-4++ classes, 
respectively (Table 5). These unknowns accounted for 75, 70, 
and 69% of the total metabolites changing significantly in 
the WT++, atferl-3-4, and atferl-3-4++ classes, respectively, 
when compared with the WT. 

The clustering of metabolites was studied in each sam- 
ple class using PLS analysis, including both the identified 
and unknown metabolites (Supplementary Fig. S2B at JXB 
online). The WT, WT++, atferl-3-4 , and atferl-3-4 ++ 
classes separated in clusters. The first vector (vl) explained 
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16% of the variability, with the second vector (v2) explain- 
ing 9% of the variability. The separation between clusters was 
associated with those metabolites with a large contribution 
(X-weight) to vl (Supplementary Table S8). Approximately 
70% of the metabolites with a high contribution to cluster 
separation were unknowns (Supplementary Table S8). The 
separation of the classes agreed with that found in the prot- 
eomic study, with atferl-3-4 and WT++ being closer together 
than WT and atferl-3-4++. 

A separate comparison was made between the triple fer- 
ritin mutant plants grown in Fe excess (atferl-3-4++) with 
those of the same mutant grown in control conditions (atferl- 
3-4 ). Fifteen identified metabolites and 21 unknowns showed 
significant changes between these two classes (Table 5). 
The largest changes in abundance were increases for urea, 
suberyl glycine, and proline (5-, 3-, and 2-fold, respectively). 
The unknowns accounted for 58% of the total metabolite 
changes. The clustering of metabolites was studied in each 
sample class using PLS analysis, including both the identified 
and unknown metabolites. The atferl-3-4 and atferl-3-4++ 
classes were well separated in clusters, with the first vector 
(vl) explaining 28% of the variability, and the second vector 
(v2) explaining 14% of the variability (Supplementary Fig. 
S2C at JXB online). The metabolites with a large contribu- 
tion (X-weight) to vl are shown in Supplementary Table S9. 
Approximately 62% of the metabolites with a high contribu- 
tion to cluster separation were unknowns (Supplementary 
Table S9). 

A separate comparison was made between atferl-3- 
4++ versus WT++. Four identified metabolites and seven 
unknowns showed significant changes between these two 
classes (Table 5). The largest changes in abundance were 
increases for urea, 2-hydroxyglutaric acid, and threonic 
acid (5-, 3-, and 2-fold, respectively), and decreases in thre- 
ose (2-fold). The unknowns accounted for 64% of the total 
metabolite changes. 

The information on the changes in metabolite abundances 
occurring in the different treatments was also clustered using 
Venn diagrams similar to those obtained for transcriptomic 
and proteomic data (Fig. 4). The category containing over- 
laps in the comparisons of atferl-3-4 versus WT, atferl-3- 
4++ versus WT++, and atferl-3-4++ versus WT has four 
metabolites: 2-hydroxyglutaric acid, threonic acid, and two 
unknowns. The second group of interest (only in atferl-3- 
4++ versus WT++) contained four metabolites: threose and 
three unknowns. The third group of interest (only in atferl-3- 
4++ versus WT) contained 23 metabolites: gluconic acid lac- 
tone, tyrosine, 3,6-anhydrogalactose, iV-acetyl-D-hexosamine, 
salicin, ascorbic acid, and 17 unknowns. 

Comparison of the flower hormonal profiles of Col-0 
and atfeii -3-4 Arabidopsis plants grown under iron- 
sufficient or excess conditions 

The atferl-3-4 mutant plants have their flower development 
severely impaired under Fe excess conditions (Ravet et al, 
2009a), suggesting some modifications in hormone con- 
centration. The hormonal changes observed in this study 



between the two genotypes under two Fe conditions were 
restricted to the cytokinin concentrations, with the changes 
being especially intense in the mutant plants grown under Fe 
excess conditions. The IAA concentration in flower extracts 
was not significantly different in the atferl-3-4 and WT plants 
grown at both Fe levels. A similar conclusion can be drawn 
for the ABA concentrations in flower extracts (Table 6). In 
contrast, the concentrations of some cytokinins, including Z, 
tZR, cZR, and DHZR, showed significant changes (Table 6), 
whereas those of iPR did not change (the iP concentrations 
were below the limit of quantification; data not shown). The 
Z concentration in flowers was significantly higher in atferl- 
3-4++ plants than in the other three plant classes, which had 
similar Z concentrations. On the other hand, similar tZR con- 
centrations were observed in flowers from WT and mutant 
plants when grown under control conditions. However, the Fe 
excess treatment led to a significant and similar increase in the 
tZR concentration in flowers of both genotypes (Table 6). The 
cZR concentration was similar in flowers of WT, atferl-3-4, 
and WT++ plants, but it was significantly lower in the flowers 
of atferl-3-4++ plants. The DHZR concentration was lower 
in flowers of atferl-3-4 plants than in the WT plants, under 
both control conditions and Fe excess (Table 6). In flowers 
from WT plants, the DHZR concentration was lower with Fe 
excess, but the decrease was not statistically significant. 

Discussion 

The role of plastids in Fe homeostasis is an important param- 
eter for normal flower development and for plant fertility. 
To document this point further, a global analysis (transcrip- 
tome, proteome, metabolome, and hormone concentrations) 
of Arabidopsis flowers from WT and triple atferl-3-4 ferritin 
mutant plants grown under control or Fe excess conditions 
was carried out. This approach proved to be valid based on 
the expression of the ferritin genes observed on the Affymetrix 
microarrays (Fig. 1C). The very low level of expression of 
the AtFerl and AtFer4 genes in the flowers of the atferl-3-4 
mutant plants confirms that these two genes were knocked 
out. On the other hand, only a mild decrease of the expres- 
sion of the AtFer3 gene in the atferl-3-4 mutant was observed 
when compared with WT plants under control conditions. 
This decrease in AtFer3 transcript abundance becomes sig- 
nificant when the atferl-3-4 mutant plants are challenged with 
an excess of Fe. This result confirms that the AtFer3 gene was 
knocked down, because a T-DNA insertion within its pro- 
moter region abolishes the regulation of its expression by Fe 
excess, but does not prevent its basal level of expression (Ravet 
et al, 2009a). The robustness of the microarray analysis was 
further assessed by re-analysing 12 representative genes by 
qRT-PCR (Supplementary Table SI at JXB online). Only one 
gene for one comparison showed deviating behaviour between 
qRT-PCR and microarray analysis. 

Changes in gene expression at the transcript level because 
of mutations within ferritin genes and/or Fe excess treatment 
were low, with a total of 276 genes being significantly regu- 
lated (Fig. 1 A, B; Table S5 at JXB online). Most fold changes 
observed in transcript abundance of these genes were not very 



Arabidopsis flower-omics of a ferritin mutant | 2681 

Table 5. Main changes in Arabidopsis flower metabolite levels. The corresponding response ratios, defined as the level in the WT++, 
atfer1-3-4, and atferl-3-4++ classes divided by the level in the WT, are shown. When the response ratio was <1 , the inverse was taken 
and the sign changed. Only metabolites showing changes in response ratios >1 .5-fold and statistically significant (at P < 0.05) are 
included in the table. 



WT++ versus WT 



atfer1-3-4 
versus WT 



atfer1-3-4++ 
versus WT 



atfer1-3-4++ 
versus atfer1-3-4 



atfer1-3-4++ 
versus WT++ 



Known metabolites 
Urea 

Suberyl glycine 
2-Hydroxyglutaric acid 
Proline 

Threonio acid 

Gluconic acid lactone 

Homoglutamine 

Phenylalanine 

Tyrosine 

Palmitic acid 

Methylhexadecanoic acid 

3,6-Anhydrogalactose 

Heptadecanoic acid ester 

Stearic acid 

Levoglucosan 

Glycerol 

Glucose 

W-Acetyl-D-hexosamine 

Pelargonic acid 

Maltose 

Malonic acid 

Threose 

Salicin 

Ascorbic acid 

Unknowns 

446831 

225867 

310413 

268365 

443451 

235436 

208850 

407649 

446856 

21 2235 

446777 

446849 

31 0006 

404568 

211956 

201005 

211932 

446884 

31 0006 

21 5062 

211896 

223865 

214528 

202834 

446842 

223629 

204448 

268353 

231248 



2.0 
1.8 
1.5 



-4.8 



12.5 



2.4 



3.8 
3.5 



1.6 



1.5 
1.5 



9.8 

2.5 

2.0 
2.0 
1.5 

-1.5 
-1.6 



-1.7 



-4.1 
-5.0 

20.6 

3.8 

3.0 

2.5 
2.3 
2.0 



1.8 

1.8 
1.7 

1.6 
1.6 

1.5 
1.5 



-1.5 



4.9 
2.5 

2.1 



1.5 

-1.5 
-1.5 

-1.6 
-1.6 
-1.6 
-1.6 
-1.6 
-1.7 
-1.7 
-1.8 
-1.8 



19.5 

2.7 
2.4 
1.9 



1.6 



-1.5 
-1.5 
-1.5 
-1.5 
-1.6 
-1.6 
-1.6 
-1.6 



4.7 
2.7 



-2.0 



3.0 
2.3 



1.5 
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Table 5. Continued. 





WT++ versus WT 


atfer1-3-4 
versus WT 


atfer1-3-4++ 
versus WT 


atfer1-3-4++ atfer1-3-4++ 
versus atferl -3-4 versus WT++ 


427262 








-1.6 


299396 








-1.6 


213194 






-1.7 




447146 






-1.8 


-1.9 


212251 






-1.8 




446633 








-1.8 


227724 








-1.8 


21 7809 






-1.9 




310785 






-2.0 




213182 








-2.2 


213198 






-2.5 


-1.9 


318412 








-2.5 


213182 








-2.5 


305055 






-2.6 




301584 








-3.3 


222093 




-3.7 






213230 


-3.8 


-3.1 


-5.0 





high, ranging from 1.4- to 3.5-fold, with only a few strong 
alterations higher than 7-fold. Such a moderate quantita- 
tive impact of mutations and/or treatment was also observed 
for proteins (63 proteins with significant changes in relative 
abundance; Table 4) and metabolites (70 metabolites with sig- 
nificant changes in relative abundance; Table 5), with also a 
few strong alterations. Of the 63 protein spots identified by 



two-dimensional IEF in Table 4, 61 corresponding genes were 
represented on the ATH1 chip used for microarray analy- 
sis. When checking transcript signal values for these genes, 
only two of them [Atlg06260, a cystein proteinase-like gene 
(CysP), and At3g25050, a xyloglucan endotransglucosylase 
(XTH3)] showed differential expression, which was signifi- 
cant for the atferl-3-4 versus WT comparison and elevated 



(A) atfer1-3-4 vs. WT atfer1-3-4++ vs. WT++ (B) atfer1-3-4++ vs. atfer1-3-4 WT++ vs. WT 
20: {} 2^18 12:^11^1 22: A 2 0 f 2 11: {} 0 ^11 

10:^8^2 11: 0 8 f 3 36: {) 9 | 2 7 4:^2^2 





atfer1-3-4++ vs. WT 
21: {} 9 J12 

32: A 18 J 14 

Fig. 4. Venn diagrams indicating major changes in Arabidopsis flower proteins and metabolites. Overlapping regulation is shown by 
Venn diagrams. Numbers indicate protein species (bold) and metabolites (regular) with significant and >2-fold changes in response 
ratios (i.e. abundance); up- and down-accumulation is indicated by orange and green arrows, respectively (for fold changes see Table 4). 
(A) Species showing differences in relative abundance in flowers eft atferl -3-4 versus WT, atferl -3-4++ versus WT++, and atfer1-3-4++ 
versus WT. Groups of interest discussed in the text are highlighted by a grey background. (B) Species showing differences in relative 
abundance in flowers of in atferl -3-4++ versus atferl -3-4 and WT++ versus WT. 
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Table 6. Hormone concentrations (IAA, ABA, Z, tZR, cZR, DHZR, and iPR) in flower extracts from Arabidopsis WT and atfer? -3-4 
plants grown in normal and Fe excess conditions. The results are presented as the means (for n=3) ±SD. 



Hormone LQ 
(pmol g _1 FW) 






WT (pmol g 1 FW) 


atferl -3-4 (pmol g" 1 FW) WT++ (pmol 
g" 1 FW) 


atfer1-3-4 ++ 
(pmol g _1 FW) 


IAA 


3 




77.76 + 3.11 


79.85 + 12.84 


81.68 + 14.62 


76.18 + 6.66 


ABA 


1.5 




165.10 + 16.16 


170.99 + 49.83 


147.07 + 29.09 


137.39 + 39.24 


Z 


0.23 




0.55 + 0.06 


0.57 + 0.05 


0.50 + 0.01 


0.68 + 0.08 


tZR 


0.09 




2.69 + 0.12 


2.74 + 0.19 


3.20 + 0.27 


3.13 + 0.29 


cZR 


0.09 




1.44 + 0.07 


1.50 + 0.07 


1.46 + 0.29 


1.26 + 0.12 


DHZR 


0.28 




0.45 + 0.08 


0.29 + 0.11 


0.33 + 0.15 


BLQ 


iPR 


0.04 




1.06 + 0.26 


1.13 + 0.34 


1.19 + 0.08 


1.29 + 0.04 


ANOVA source 


IAA 


ABA 


Z 


tZR cZR 


DHZR 


iPR 


Plant class 


NS 


NS 




NS NS 




NS 


Fe condition 


NS 


NS 


NS 


NS 


NS 


NS 


Plant classxFe condition 


NS 


NS 








NS 



WT, Arabidopsis Col-0 grown in normal conditions; atfer1-3-4, Arabidopsis mutant grown in normal conditions; WT++, Arabidopsis Col-0 
grown in Fe excess; atfer1-3-4++, Arabidopsis mutant grown in Fe excess. 

IAA, indole-3-acetic acid; ABA, abscisic acid; Z, frans-zeatin; tZR, frans-zeatin riboside; cZR, c/s-zeatin riboside; DHZR, dihydrozeatin 
riboside; iPR, A/ 6 -isopentenyladenosine; LQ, limit of quantification; BLQ, below the limit of quantification. 

Statistical significance of the interaction among the different parameters: one-way ANOVA, LSD Fisher test (^significant differences P < 0.05). 



for atferl-3-4++ versus WT++. To verify transcriptomic 
data, mRNA contents of CysP and XTH3 were additionally 
analysed by qRT-PCR as described in Supplementary Table 
SI. Significance analysis by a Student's ?-test revealed the 
same regulation as for DNA microarrays (data not shown). 
Overall, the multivariate statistical analysis (Supplementary 
Fig. S2) from proteomic and metabolomic data suggests that 
the mutant in control conditions clusters closely to WT plants 
grown in excess Fe, in agreement with the lower ferritin and 
higher Fe concentrations in flowers reported for this mutant 
(Ravet et al, 2009a). 

Functional categories consistently affected at the three 
levels of gene expression 

In a recent study, Amiour et al, (2012) examined the impact 
of N deficiency in maize at the transcriptomic, proteomic, and 
metabolomic levels. These authors observed that integration 
of the three 'omics' studies was not straightforward, since dif- 
ferent levels of regulation seem to occur in a stepwise manner 
from gene expression to metabolite accumulation. The same 
conclusion can be reached from the present study. Although 
no direct overlaps between transcripts, proteins, and metab- 
olites were detected, some major modifications in specific 
functional categories were, however, consistently observed at 
these three levels of expression. The strongest impact concerns 
redox reactions and oxidative stress (i) with 70 transcripts of 
genes belonging to this functional category among the 276 
affected (Fig. 1A, B; Supplementary Table S5 at JXB online); 
(ii) with changes in abundance of four proteins directly related 
to stress [monothiol glutaredoxin-S15, PR5, glutathione 
^-transferase (GST6), and peroxidase 17], and several oxidore- 
ductases (Table 4); and (iii) with a decrease in ascorbic acid, 
which appears more dependent upon Fe excess treatment than 
upon the genotype considered (Table 5). These observations 
are consistent with the various flower phenotypes reported for 



atferl -3-4 plants (Ravet et al, 2009«), which indicated that an 
Fe-mediated oxidative stress was enhanced when the ferritin 
Fe buffering capacity in chloroplasts was disturbed, ultimately 
leading to flower damage and fertility loss (Ravet et al. , 2009a; 
Briat et al, 20106). The second functional category being con- 
sistently affected at the three different levels of gene expression 
is related to amino acid and protein metabolism. The focus of 
regulated proteins in this category (Table 4) is well reflected 
by transcript regulation under Fe excess {atferl -3 -4++ versus 
WT++; Table 2). Proteolysis seems to play an important role, 
since genes coding for peptidase proteins are regulated (Tables 
1 , 2) and a cysteine proteinase-like protein (CysP) is synthe- 
sized de novo (100-fold increase) in the atferl-3-4 flowers from 
plants treated with excess Fe (Table 4). Indeed, analysis of the 
expression of the gene encoding CysP by qRT-PCR revealed 
increased transcript ratios in atferl-3-4 versus WT and atferl- 
3-4++ versus WT++ (data not shown). Further, a new pro- 
tein spot of a proteasome subunit, appearing in atferl-3-4++ 
versus atferl-3-4 profiles (Table 4), points to protein degrada- 
tion upon Fe stress in ferritin mutants. Because proteolysis- 
associated genes also represent targets for gene regulation 
in knockout and overexpressing mutants of the chloroplast 
Fe-transport protein PIC1 (Duy et al, 2007, 2009), this func- 
tional category is closely linked to plastid Fe homeostasis and 
transport. At the metabolite level, the accumulation of urea 
(Table 5) is the ultimate marker of this activated proteolysis. 
Also consistent with changes in protein metabolism is the 
alteration in the pools of some amino acids such as tyrosine, 
phenylalanine, and proline in atferl-3-4 mutant plants treated 
or not with excess Fe (Table 5). 

Influence of the ferritin mutation background on the 
flower-omes regardless of external iron feeding 

The Venn diagrams (see shading in Figs 1 A and 4A) show that 
the expression of 54 genes (Table 1) was changed, whereas no 
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proteins and only four metabolites were affected exclusively 
by the effect of the mutations. Among the genes affected, the 
stress and redox categories accounted for 28% of the changes. 
They also contained those genes showing the largest differ- 
ences in expression. The majority of them showed up-regu- 
lation, some of them encoding Fe-haem oxidoreductases 
from the CYP family (cytochrome P450 superfamily), FAD 
berberine family, and two cell wall peroxidases. The substrate 
of CYP includes lipids, hormones, and xenobiotic substances, 
while the FAD-binding berberine family are, in general, 
membrane-associated oxidoreductases involved in electron 
transport, with both families usually being expressed dur- 
ing flower anthesis (TAIR annotation). The absence of haem 
proteins in the 2-DE study is probably related to the fact that 
most of them are membrane-located or cell wall-associated 
proteins, and therefore difficult to detect by this technique. 
Transcriptomic data also indicate that the impact of the muta- 
tions in the redox homeostasis processes occurs mainly at the 
membrane level. Regarding the metabolites found in this over- 
lap, the increase in threonic acid, a product of ascorbic acid 
catabolism (Helsper and Loewus, 1982), would also suggest an 
increased redox stress. The increase in 2-hydroxyglutaric acid 
is also in agreement with the alterations caused by the mutant 
background, since it increases during senescence, when prote- 
olysis and P-oxidation occur (Engqvist et al, 2009, 2011). 

Another category strongly affected by the mutations is 
hormone metabolism (Table 1), which shows a general up- 
regulation of genes responsive to ABA (At5g42560), ethylene 
(Atlg36060, At3g50260), and auxin (At3gl7600, Atlg05680). 
The latter gene encodes a UDP-glucosyltransferase that acts 
on indole-3-butyric acid (therefore affecting auxin homeosta- 
sis), whose transcript and protein levels are strongly induced 
by H 2 0 2 , and which may allow integration of ROS and 
auxin signalling. In consonance with these results, it is worth 
mentioning that an additional target for gene regulation 
contained two transcription factors, At5g46350 (WRKY8) 
and At2g43000 (ANAC042, see Table 1), which are usually 
induced by the presence of H 2 0 2 , supporting the close connec- 
tion between redox homeostasis and the mutant phenotype. 

Impact of ferritin mutations on the flower-omes under 
iron excess in comparison with the wild type 

The Venn diagrams also revealed 74 genes, eight proteins, 
and four metabolites differentially regulated only in atferl-3- 
4++ versus WT++, which mainly reflect the consequences of 
excess Fe on the mutant (see shadings in Figs 1A, 4A; Tables 
2, 4, 5). Among these genes, the specific increase in the FR03 
transcripts is surprising since it is described to be up-regulated 
in response to Fe deficiency. It may indicate amplified Fe(III) 
reduction activity and/or more transport and distribution of 
Fe. Targets for gene regulation in additition to the aforemen- 
tioned proteolysis function also included the stress and redox 
and hormone metabolism categories. However, the effects of 
Fe excess (Table 2) were, in general, opposite to those exerted 
exclusively by the mutation regardless of the Fe status (Table 1). 
More specifically, FAD oxidoreductases and peroxidases were 
down-regulated, whereas CYP genes were still up-regulated. 



This different response of the Fe-haem proteins may point 
towards their different roles in the Fe and/or redox homeo- 
stasis mechanisms. In agreement with the transcriptomic data, 
two oxidoreductase proteins (spots 20 and 51) increased in 
abundance in the same comparison (Table 4). General stress 
responses represented by pathogenesis-related and disease 
resistance genes and by GST6 in the proteomic study (spot 
40, Table 4) were also up-regulated. This up-regulation of 
stress responses, not observed in the WT, may indicate that 
atferl-3-4++ flowers are facing more intense stress conditions 
than those of the WT++. In this comparison (Table 2), genes 
encoding hormone-responsive genes, similar to those up-regu- 
lated as a consequence of the mutations (Table 1), were down- 
regulated and included two auxin- (Atlgl9830, At2g46530) 
and two ethylene- (Atlg63030, At5g52020) responsive genes. 
In addition, three of the proteins showing an increase in rela- 
tive abundance in this group (spots 40, 48, and 51) are auxin 
regulated (UniProt DB annotation). 

Adaptation of flower-omes representing an Fe effect 
on the ferritin mutant 

Comparisons of atferl-3-4++ versus WT showed changes 
in 43 transcripts, 14 proteins, and 23 metabolites (shading in 
Figs 1A, 4A, Tables 3, 4, 5), reflecting the Fe effect on the 
mutant. Besides the already described target categories for 
regulation, chloroplast-encoded genes are strongly regulated 
(Table 3; Supplementary S5 at JXB online) at the transcript 
level, especially in the atferl-3-4++ versus WT comparison. 
This is consistent with three photosynthesis proteins affected 
(Table 4). It may indicate that pronounced imbalance in chlo- 
roplast Fe homeostasis by loss of the plastid Fe buffer fer- 
ritin, in parallel with excess external Fe supply, leads to the 
reprogramming of some chloroplast-intrinsic processes. Also, 
the strong decrease of salicin and ascorbic acid (Table 5), 
which are synthesized (salicin) and function as an antioxidant 
(ascorbic acid) in chloroplasts, might reflect disturbances of 
the plastid metabolism. 

Cell wall and membrane biology modifications in 
response to ferritin mutations and/or Fe excess 

Target categories for gene regulation by Fe excess on the 
mutant versus WT also included cell wall modification and 
transport, which accounted for 12% and 8% of the changes, 
respectively (Table 2). At the transcript level, a down-regu- 
lation was observed in genes encoding two pectin esterases 
and several structural components of the cell wall, including 
two expansins (EXPB4 and EXPA22), FLA12 (an arabinoga- 
lactan protein), and GRP5 (glycine-rich protein 5), whereas 
an up-regulation was observed in two extensins (Atlg26240, 
At3g28550). At the protein level (Table 4), there was an 
increase in xyloglucan endotransglucosylase (spot 48), a cell 
wall-loosening enzyme thought to promote cell expansion 
(Fry et al., 1992). When re-analysed by qRT-PCR, increased 
transcript ratios were found for the xyloglucan endotransglu- 
cosylase (XTH3) in atferl-3-4 versus WT and atferl-3-4++ 
versus WT++ comparisons (data not shown). These changes 
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in cell wall structure may be behind the impaired fertility 
observed in the mutant, especially when grown with Fe excess. 
Furthermore, several of these genes are auxin responsive, and 
this might be directly related to the complex hormone bal- 
ance commented on above. 

Regarding transporter genes, 12 of them showed changes 
in regulation because of the atferl-3-4 mutations and/or Fe 
excess. Nitrate (NTR1.1) and phosphate (PHT4.3) transport- 
ers were up-regulated in the mutant background (Table 1), 
whereas an array of transporters were changed by Fe excess 
(Tables 2, 3), including several members of the ABC fam- 
ily (ALS3, ABCG10, and ABC13), a flavonoid antiporter 
(TT12), the amino acid (CAT6), malate (ALMT1), uracil 
(UPS2), and sugar (ZIFL2) transporters, and two transport- 
ers of unknown specificity (DTX3 and SEC14), As expected 
from their membrane location, no transporters were found in 
the proteome analysis, whereas changes in the metabolome 
indeed reflected the changes in transporters described above. 
In particular, the increases in the relative abundance of amino 
acids (proline and phenylalanine) and the decreases in the 
relative abundances of sugars (maltose and glucose) would 
be in agreement with the increases observed in CAT6 (amino 
acid transporter) and decreases in ZIFL2 (sugar transporter) 
transcript abundances in the mutant. Further, the focus on 
regulation of transport-associated genes and the high overlap 



(about one-third, see Table 2) with genes also differentially 
expressed in PIClox lines (Duy et al, 2011) indicate a close 
correlation between disturbance of Fe homeostasis and mem- 
brane transport. A strong deregulation of several transport- 
ers had already been reported in atferl-3-4 plants (Ravet 
et al, 2009a). The present results are not fully consistent 
with previous results, with regard to the Fe(II) transporter 
IRT1 (Vert et al, 2001), and the chloroplast iron permease 
PIC1 (Duy et al, 2007). In the present microarray study, also 
confirmed by qPCR (data not shown), IRT1 and PIC1 did 
not show changes in expression between WT and atferl-3- 
4 plants, whereas IRT1 was found to be strongly down-reg- 
ulated (Ravet et al, 2009a) and PIC1 overexpressed (Duy 
et al, 2007) in the atferl-3-4 mutant flowers. The reason(s) 
for such discrepancies is so far unknown. 

Data discussed above indicate that mutations in ferritin 
genes exert an effect on redox enzymes, especially in those 
located in membranes, and that both ferritin mutations and 
Fe excess affect membrane transporters. These effects on 
membrane biology also correlate with changes observed in 
the lipid metabolism at the three omic levels. At the tran- 
scriptome level, the mutant background caused a decrease 
in a GDSL-motif lipase and in a lipid transferase (Table 1), 
which is consistent with the decreases measured in the rela- 
tive amount of two GDSL lipases (spots 36 and 37, Table 4) 
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Fig. 5. Schematic representation of the biological functions bridging chloroplast Fe homeostasis to flower development and fertility. 
Disturbing Fe chloroplast homeostasis by altering the ferritin content promotes an oxidative stress and alters sugar, amino acid, and 
protein metabolism. This is consistent with the observation that chloroplast functions are perturbed, and that among transporters 
whose expression is altered in the atferl-3-4 ferritin mutant some of them are sugar and amino acid transporters. The membrane (lipid 
metabolism) and cell wall structure are also modified, as well as the hormonal balance, including cytokinins and perhaps auxins. This 
synthetic representation is based on the comparison (between Col-0 and the atferl-3-4 ferritin mutant plants grown under standard or 
excess Fe conditions) of the changes in the transcriptome, proteome, and metabolome occurring in the various functional categories 
defined in the text. 
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at the proteomic level, and at the metabolomic level with the 
changes in the pool of fatty acids (pelargonic, heptadecanoic, 
palmitic, and methyhexadecanoic acids), especially in atferl- 
3-4++ versus atferl-3-4 (Table 5). 

Hormonal imbalance in response to ferritin mutations 
and/or Fe excess 

Among the various hormone concentrations measured in 
this study, only the concentration of active Z-based cyto- 
kinins was increased under Fe excess conditions in flowers of 
both genotypes (Table 6). This effect was more pronounced 
in flowers of atferl-3-4 mutant plants. Accordingly, in atferl- 
3-4++ versus WT (Table 3), a down-regulation of the mito- 
chondrial glutamate dehydrogenase (GDH3, At3g03910) was 
observed. This enzyme has been shown to be up-regulated 
in response to cytokinin and probably involved in control of 
nitrogen metabolism (Igarashi et al, 2009), which also seems 
changed in the mutant, as indicated by the increase in urea 
in the atferl-3-4++ versus WT and atferl-3-4 comparisons. 
A link between Fe homeostasis and cytokinins has already 
been reported for the regulation of the root Fe uptake system 
(Seguela et al, 2008). It is well known that the main active 
forms of cytokinins tended to concentrate in those plant 
organs that have high metabolic activity (Marschner, 1995; 
Sakakibara, 2006). The presence of cytokinins in these meta- 
bolic sites has a sink effect, attracting those metabolites and 
nutrients that are relevant to the development of intense met- 
abolic activity (Marschner, 1995; Sakakibara et al, 2006). In 
the context of the present experiments, the increase in concen- 
tration of Z-based active cytokinins in flowers suggests that Fe 
excess might increase the need for metabolites and nutrients 
in order to cope with the stress situation. Further, the up-reg- 
ulation observed in nitrate and phosphate transporters in the 
mutant background would be in agreement with this hypoth- 
esis (Table 1). This situation seems to be more relevant in the 
case of atferl-3-4 mutant plants. This is supported by the fact 
that the increase in tZR is associated with a decrease in cZR, 
thus suggesting the conversion of cZR into tZR through the 
induction of cis-trans isomerase activity (Sakakibara et al. , 
2006). Other hormones have been described to play a role in 
Fe metabolism, in particular auxin and ethylene (Bacaicoa 
et al, 2011; Lingam et al, 2011; Romera et al, 2011; Giehl 
et al, 2012), and also in flower functionality, as it is the case 
for gibberellins and IAA (Davies, 2004). Although a large 
number of genes and proteins regulated by these hormones 
had changed, the present study did not show any differences 
in the concentration of these plant hormones in flowers of 
both WT plants and atferl-3-4 mutant plants. However, since 
the analysis quantified the total hormone concentration, 
the possibility that conjugated forms or cellular distribution 
could be altered by Fe treatment cannot be excluded. 

Influence of the Fe excess treatment on the flower- 
omes regardless of the genotype 

Finally, six genes, one protein, and one metabolite respond- 
ing to excess Fe irrespective of the genotype (common in 



the WT++ versus WT and atferl-3-4 ++ versus atferl-3-4 
comparisons; Figs IB, 4B) could be pinpointed. Among 
these genes, two encoded heat shock proteins and one a 
receptor kinase, and this is consistent with a stress situa- 
tion. Unfortunately, two of these genes (Atlg20070 and 
Atlg76240, Table 3) and the protein (spot 56, Table 4) were of 
unknown function, and the metabolite (446831, Tables 4, 5) 
was also an unknown, but interestingly it was the one show- 
ing the largest change in the metabolomic study. These eight 
Fe-responsive elements are attractive targets for future stud- 
ies dealing with the plant responses to Fe toxicity. 

Conclusions 

In conclusion, disturbing Fe homeostasis by mutations in fer- 
ritin genes and/or by Fe excess treatment leads to alteration of 
flower biology. This is reflected by changes in gene expression 
at the transcript and protein levels, as well as by modifica- 
tions in various metabolite pools. Functional categories con- 
sistently affected at those three levels of gene expression were 
oxidative stress and protein metabolism. The mutant back- 
ground caused alterations in Fe-haem redox proteins located 
in membranes, as denoted by changes in their transcript 
abundances and in hormone-responsive proteins (Fig. 5). 
The specific effect of excess Fe in the mutant caused further 
changes in these categories, suggesting that the mutant is fac- 
ing a more intense Fe/redox stress. Changes in the membrane 
transporters and lipid metabolism (Fig. 5) also indicate that 
the impact of the mutation and/or Fe excess are very intense 
at the membrane level. In spite of the large number of genes 
and proteins responsive to hormones found to be regulated 
in this study, changes in the hormonal balance are restricted 
to cytokinin concentrations, being especially obvious in the 
mutant plants grown under Fe excess conditions. This reflects 
the complexity of the hormonal balance and the wide range 
of targets for hormone action in flower metabolism. 

Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Protein species identified in extracts of 
Arabidopsis flowers. 

Figure S2. Partial least square (PLS) analysis of Arabidopsis 
flower protein species and metabolites. 

Table SI. Verification of DNA microarray data by quanti- 
tative RT PCR. 

Table S2. Overview of the metabolite reporting list. 

Table S3. Optimal parameters for the detection of IAA and 
ABA used in the HPLC-ESI-MS/MS analysis of Arabidopsis 
plants. 

Table S4. Optimal parameters for the multiple reaction 
monitoring detection of CKs used in the HPLC-ESI-MS/MS 
analysis of Arabidopsis plants. 

Table S5. Genes significantly regulated in all five compari- 
sons of atferl-3-4 mutants versus Col-0 wild type (WT) under 
Fe-sufficient and Fe -excess (++) conditions. 

Table S6. Summary of the 2-DE protein profiling study 
results. 
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Table S7. Complete list of identified and unknown 
metabolites. 

Table S8. Metabolite X-weights in the partial least square 
(PLS) analysis of Arabidopsis flowers shown in Supplementary 
Fig. S2. 

Table S9. Metabolite X-weights in the partial least square 
(PLS) analysis of atfer 1-3-4 flowers shown in Supplementary 
Fig. S2. 
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